Studying Exotic Hadrons in Heavy Ion Collisions 
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We investigate the possibilities of using measurements in present and future experiments on heavy 
ion collisions to answer some longstanding problems in hadronic physics, námely identifying hadronic 
molecular states and exotic hadrons with multiquark components. The yields of a selected set of 
exotic hadron candidates in relativistic heavy ion collisions are discussed in the coalescence model 
in comparison with the statistical model. We find that the yield of a hadron is typically an order of 
magnitude smaller when it is a compact multiquark statě, compared to that of an excited hadronic 
statě with normál quark numbers. We also find that some loosely bound hadronic molecules are 
formed more abundantly than the statistical model prediction by a factor of two or more. Moreover, 
due to the significant numbers of charm and bottom quarks produced at RHIC and even larger 
numbers expected at LHC, some of the proposed heavy exotic hadrons could be produced with 
sufficient abundance for detection, making it possible to study these new exotic hadrons in heavy 
ion collisions. 

PACS numbers: 14.40.Gx,11.55.Hx,24.85.+p 



I. INTRODUCTION 

Experiments at the Relativistic Heavy Ion Collidcr 
(RHIC) during the past decadc háve shown that the prop- 
erties of the quark-gluon plasma (QGP) formed in heavy 
ion collisions are far more intriguing than originally con- 
ceived |l| . Instead of weakly interacting, the quark-gluon 
plasma was found to be a strongly coupled systém with so 
small a shear viscosity that it behaves like an ideál fluid. 
The study of the QGP is expected to remain an active 
field of research in the future because of the proposed 
upgrade of RHIC and new experimental possibilities at 
the Large Hadron Collidcr (LHC). The physics of QGP 
is related to a wide range of other fields, such as the early 
universe, nonequilibrium statistical physics, string theory 
by AdS/CFT correspondence, and so on. 

During the samé time, there háve also been exciting de- 
velopments in the spectroscopy of heavy hadrons, start- 
ing with the discovery of the £> sJ (2317) by the BABAR 
Collaboration Q, whose mass could not be explained by 
the simple quark model, and the X(3872) by the Belle 
Collaboration [Š| , whose mass and decay channel strongly 
support a nontrivial fraction of DD* and DD* compo- 
nents in its wave function. The Belle collaboration also 
reported the fmding of the Z + (4430) in the vb'-K + spec- 
trum J4J. Because the Z + (4430) is a charged statě, it 



cannot be a simple cč statě. If confirmed, this would be 
the hrst evidence for the existence of an exotic hadron 
that is composed of two quarks and two antiquarks like 
cčuď @ . 

The question of whether multiquark hadrons exist is 
an old problém in the light hadron sector that began 
with attempts to understand the inverted mass spec- 
trum of the scalar nonet (on(980), /o (980), and so on) 
in the tetraquark picture 0, [D]- Also, the exotic H 
dibaryon was proposed on the basis of the color-spin 
interaction [8j, and it has been sought for in various 
experiments for a long time without success. On the 
other hand, the A(1405) baryon resonance was considered 
as a KN quasi-bound statě in the KN — 7rS coupled- 
channel analysis [9J]. With further development of the 
coupled-channel approach [lOJ . it has been realized that 
the A(1405) is the most obvious and uncontroversial can- 
didatc for a hadronic moleculc, whose wave function is 
composed dominantly of a KN bound statě mixed with 
a small 7rE resonant statě 11|. If such a configura- 
tion for the A(1405) is confirmed in experimental mea- 
surements, it will be the first evidence for a molecular 
hadronic statě. Many new multiquark states such as 
KKN, KNN, (ížíž)o, and etc. háve also been predicted. 

Whilc the Z + (4430) could be a first explicitly exotic 
hadron found to datc, it will be a milestone in hadronic 



spectroscopy if other flavor cxotic hadrons are found, 
such as the controversial pentaquark <d + (ududš) first re- 
ported in experiments on photonuclear reactions J12J . 
Also, several other flavor exotic molecular and compact 
multiquark hadrons were previously predicted, based, re- 
spectively, on the meson-exchange and color-spin inter- 
actions. More recently, a sim ple diquark model based on 
the color-spin interaction [13l Il4| was shown to naturally 
cxplain the likely existence of flavor exotic multiquark 
hadrons consisting of heavy spectator quarks, such as the 
T^ c {udčč), T° b (udčb), and Q cs (udusč) that were predicted 
before, and the newly predicted dibaryon H^ + (udusuc) . 
Furthermore, the stable molecular bound states ĎN and 
DNN in the charmed sector and the BN and BNN in 
the bottom sector háve been predicted to exist as a result 
of the long-range pion-exchangc potential [15|, [l6| . 

To gain insights into all these proposals and ques- 
tions Il3l . If 71 Il8| . we háve proposcd in recent publica- 
tions [19L I2OJ] a new approach of studying exotic hadrons 
in heavy ion collisions at ultrarclativistic cncrgies. Thcrc 
are several merits for this approach compared to the 
search for exotic hadrons in elementary particle reactions 
that háve been pursued so far. First, an appreciable num- 
ber of heavy quarks are expected to be produced in these 
collisions, reaching as largc as 20 čc pairs per unit ra- 
pidity in Pb+Pb collisions at the Large Hadron Collidcr 
(LHC) |21| . Second, through vertex reconstruction of 
weakly decaying particles, one could substantially reduce 
the backgrounds in the dctcction, making the identifica- 
tion of weakly decaying exotics possible. Finally, because 
of the large volume of quark and hadronic matters formed 
in these collisions and the new paradigm of hadronization 
through coalescence |22h24J , various exotic hadrons could 
be formed from the recombinations of quarks. We háve 
thcrcfore investigated in Refs. [f9|, Í2Q[ the possibility of 
using mcasuremcnts in present and future experiments 
on heavy ion collisions to answer the longstanding prob- 
lems in hadronic physics of identifying and examining 
hadronic molecular states and exotic hadrons that con- 
sist of multiquarks. 

In the present páper, we extend the discussions in 
Ref. [I9( to include all exotic hadrons that háve been 
proposed so far. These hadrons are classified as exotic 
mesons, exotic baryons, and exotic dibaryons, and for 
each exotic hadron we give a brief summary of its prop- 
erties and the status on the latest researches. Moreover, 
we collect all necessary references so that the present pá- 
per can be ušed as a generál guide to the literatures on 
exotic hadrons. Using these information, we then cvalu- 
ate their yields in heavy ion collisions based on both the 
statistical and the coalescence model. In particular, we 
present a detailed deseription of the calculations carried 
out in Ref. pij with an emphasis on the expected yields 
of exotic hadrons in heavy ion collisions at LHC. 

The statistical model is based on the assumption that 
hadrons produced in relativistic heavy ion collisions are 
in thermal and chemical equilibrium at temperatures 
that are close to that for the QGP to hadronic matter 



phase transition, and it has been known to deseribe very 
well the relative yields of normál hadrons J25( ■ The co- 
alescence model, on the other hand, deseribes hadron 
production through the coalescence or recombination of 
particles J22T - |24| . The model has successfully explained 
observed enhancement in the production of midrapid- 
ity bar yon s in the intermediate transverse momentům 
region [261 |27J, the quark number scaling of the ellip- 
tic flow of identified hadrons [28|, [29| , and the yield of 
antihypertritons recently discovered in heavy ion colli- 
sions at RHIC [301 ■ Furthermore, the coalescence model 
that takés into account the internal structure of hadrons 
has been shown in Ref. [31[ to be able to deseribe the 
observed suppression of the A(I520) yield in heavy ion 
collisions at RHIC compared to the prediction of the 
statistical model (32[. As in Ref. [19j . we fix the pa- 
rameters in the coalescence model by reproducing the 
yields of normál hadrons in the statistical model [25| . 
We then apply these two models to calculate the yields 
of exotic hadrons. As reference particles for the compar- 
ison with exotic hadrons, we consider, in particular, nor- 
mál hadrons such as the strange A(1115) and A(1520), 
the charmed A c (2286), and the bottom A h (5620). Our 
results show that the yields of proposed exotic mesons, 
baryons, and dibaryons are large enough for carrying out 
rcalistic measurements. We further discuss their most 
probable weak decay channels that can be observed in 
experiments. 

Our study shows that results from the coalescence 
model are sensitive to the inner structure of hadrons, 
such as the angular momentům, quark numbers and so 
on. This is different for the statistical model, because the 
yields of the hadrons in this model are determined mainly 
by their masses. To diseriminate the different pictures for 
exotic hadrons, we can thus compare the results from the 
coalescence model with those from the statistical model. 
We find that the relative yields of light exotic hadrons 
in the coalescence model are very different from those 
in the statistical model, and this makes it possible to 
experimentally diseriminate among the different pictures 
for their structures, such as multiquarks versus hadronic 
molecules. 

This páper is organized as follows. In Sec. [n] we de- 
seribe briefly the schematic model ušed for deseribing the 
time evolution of the hot dense matter formed in rela- 
tivistic heavy ion collisions. We then cxplain the statisti- 
cal model and the coalescence model ušed in the present 
study. The properties of the exotic hadrons included in 
this work and their predicted yields are given in Secs. IIIII 
and IIV1 respectively. This is followed by discussions in 
Sec. [V] and conclusions in Sec. [VT] In Appendix A, we 
derive the Wigner function of hadrons whose structures 
are deseribed by the cř-wave and the corresponding coa- 
lescence factor. We further give in Appendix B the coa- 
lescence factors for the generál Z-wave harmonie oscillator 
wavc functions. 



II. MODELS FOR HADRON PRODUCTION 
A. Heavy ion collision dynamics 

For the dynamics of centrál relativistic heavy ion col- 
lisions, we use the schematic model of Ref. [17| based on 
isentropic boost invariant longitudinal and accelerated 
transverse expansions. In this model, both the initial 
quark-gluon plasma and hnal hadronic matter are treatcd 
as noninteracting free gas, and the transition between 
these two phases of matter is taken to be first-order. The 
time evolution of the temperature and volume of the sys- 
tém is determined by the entropy conservation. In Ta- 
ble U we tabulate the values of the critical or hadroniza- 
tion temperature Tq and volume Ve at the beginning of 
the quark-gluon to hadronic matter phase transition, the 
volume Vh at the end of the mixed phase or hadroniza- 
tion, and the kinetic or thermal freeze-out temperature 
Tf and volume Vf for both centrál Au+Au collisions at 
y/IWŇ = 200 GcV at RHIC and centrál Pb+Pb collisions 
at y/sjyN = 5 TeV at LHC. Also given in Table Uis the 
abundance of various quarks at Tc ■ 



TABLE I: Quark numbers at hadronization temperature Tc 
and volume Ve, the volume Vh at the end of hadronization, 
and the thermal freeze-out temperature Tf and volume Vf in 
centrál heavy ion collisions at RHIC and LHC. 





RHIC 


LHC 


N u = N d 


245 


662 


N s = N š 


150 


405 


N C = N Č 


3 


20 


N b = N- b 


0.02 


0.8 



Ve 1000 fm 3 2700 fm 3 

T c = T„ 175 MeV 175 MeV 

V H 1908 fm 3 5152 fm 3 

ii B 20 MeV MeV 

tis 10 MeV MeV 

Vf H322 fm 3 30569 fm 3 

T F 125 MeV 125 MeV 

As deseribed in next subsections, hadron production 
in the statistical model oceurs in the volume Vh at the 
temperature Th, which is assumed to be the samé as Tc, 
while in the quark coalescence model, the production of 
both ordinary and multiquark hadrons oceurs at the tem- 
perature Tq when the volume of the QGP is Ve- For the 
production of hadronic molecular states from the coales- 
cence of hadrons, it takés pláce, on the other hand, in 
the volume Vf at the kinetic freeze-out temperature Tp. 



heavy ion collisions. In this model, the number of pro- 
duced hadrons of a given type h is given by |25J 



Nf at =Vw — 

h H 2tt2 



p 2 dp 



o lh 



l e E h /T H ± l 



Jh9hVH 2rp Tf ( irp \ 

-m h T H K2(m h /T H ) 



2tt 2 



ir Íh9hVH 



í m h T H 
V 2?r 



3/2 



i/T H 



(1) 

(2) 
(3) 



In the above equations, gn is the dcgcncracy of the hadron 
and 7/j is the fugacity, and Vh and Th are, respectively, 
the volume and temperature of the source when statisti- 
cal production oceurs. 

Since strangeness is known to reach approximatc chem- 
ical equilibrium in heavy ion collisions at RHIC due to 
the short equilibration time in the quark-gluon plasma 
and the net strangeness of the QGP is zero, the strange 
chemical potential is small and is taken to be [i s = 10 
MeV. Its value deereases with inereasing collision energy 
and is assumed to be MeV in heavy ion collisions at 
LHC. For charm and bottom quarks, they are produced 
from initial hard scattering and their numbers are much 
larger than those expected from a chemically equilibratcd 
quark-gluon plasma. As a rcsult, we obtain the fugac- 
ity Jh > 1 for both charmed and bottom hadrons. In 
terms of the fugacities 7 C and 7;, of charm and bottom 
quarks, the fugacities of charmed and bottom hadrons 
are products of 7™ and 7™ where n and m are, respec- 
tively, the charm and bottom quark numbers in these 
hadrons. Therefore, the fugacity of hadron species h can 
be written generálky as 



lh = 7™<= + "=7™ i>+ ™ E e (AIBB+ ' Í35:)/TH 



(4) 



where B, S, n c (n s ), Ji&(řig) are the baryon number, 
strangeness, (anti-)charm quark number, and (anti- 
)bottom quark number of the hadron, respectively. 

For the charm and bottom fugacities j c and 7;,, they 
can be determined by rcquiring that the total yield of 
charmed or bottom hadrons estimated in the statistical 
model is the samé as the expected total charm (N c ) or 
bottom (Ni,) quark number from initial hard nucleon- 
nucleon scattering. Using the values N c = 3 and 2V& = 
0.02 for heavy ion collisions at RHIC, we obtain j c = 6.40 
and 7b = 2.2 x 10 6 according to the following calculations: 

N c = N D +N D * + ^ (N Ds + N Ďt ) + i (N Ac + N Kc ) 

0.33 + 0.29 0.14 + 0.11 
= 1.04+1.53+ + = 3, (5) 



B. The statistical model 

The statistical model has been shown to deseribe very 
wcll the relative yields of normál hadrons in relativistic 



N b - N B + N B , + 1 (N Bs + N Bs ) + 1 (N Ab + N Ab ) 



= 5.3 x 10" 
1.7+1.5 



1.23 x 10" 



x 10" 



9.2 + 7.3 



x 10" 



0.02. (6) 



In the above, we háve ušed the average yield of heavy 
strange and antistrange mesons as well as that of heavy 
baryons and antibryons to rcmove the effect of baryon 
and strangeness chemical potentials. A similar analysis 
for heavy ion collisions at LHC based on the charm and 
bottom quark numbers N c = 20 and Nb = 0.8 (see Ta- 
ble [T| gives the charm and bottom fugacities j c = 15.8 
and 7b = 3.3 x 10 7 . These values together with those for 
RHIC are given in Table M 



TABLE II: Fugacities for c and b quarks, and hadron numbers 
from the statistical model at thermal freeze-out temperature 
Tf and volume Vf in centrál heavy ion collisions at RHIC and 
LHC, including contributions from resonance decays shown in 
the fourth column. 



coalcscence of n constituents is given by [2í 





RHIC 


LHC 




7c 


6.40 


15.8 




1b 


2.2 x 10 6 


3.3 x 10 7 




N K 


142 


363 


K,K* 


N R 


127 


363 


K,K* 


N D =N Ď 


1.0 


6.9 




N D * = N D , 


1.5 


10 




N Dl 


0.19 


1.3 




N B =N b 


5.3 x 10" 3 


0.21 




N b * = N B , 


1.2 x 10" 2 


0.49 




N N 


62 


150 


N,A 


N e 


4.7 


13 




N n 


0.81 


2.3 




N Sc 


0.10 


0.65 





In Table Hl we also show the yield N K {N R ) of K (K) 
mesons given as a sum of the directly produced K (K) 
and those from the strong decay of K* (K*) after their 
freeze-out. Similar ly, the yield of nucleons Njy includes 
both directly produced TV and those from the strong de- 
cay of A. These rcsults are obtained by using the Fcrmi- 
Dirac and Bose-Einstein distributions in Eq. ((T|). We 
notě that approximating these distributions by the Boltz- 
mann distribution as given in Eq. ([2]) does not introducc 
a large error (at most 1 % for exotic hadrons) , while the 
non-relativistic approximation ušed in Eq. (J3)) leads to an 
error of about 30-40%. 



The coalescence model 



N, 



coal 



9h 



ní pi ■ da l d 3 p, 



xf W (xi 7 - ■■ ,X n ipi,--- ,p n ). 



(7) 



In the above equation, g^ is again the degeneracy of the 
hadron whereas gi is that of its ith constituent and &Ui 
denotes an element of a space-like hypersurface. The 
function f(xi,Pi) is the covariant phase-space distribu- 
tion function of the constituents in the emission source, 
and it is normalized to their number, i.e., 



Pí ■ d<7i 



d 3 p, 

(2tt)3S 



f(xi,pi) = N i: 



(8) 



and the function f w \x\...x n : p\...p n ) is the Wigner func- 
tion of the produced hadron and is dehned by 



/ (xi,--- ,x n :pi,--- ,Pn) 



/n 
J] d yi e iPiyi iP*(xi + yi/2, ■■■ ,x n + y n /2) 

x.i/>(xi-yi/2,--- ,x n -y n /2), (9) 



; = 1 



mfc 



in terms of its wave function il){x\, • • ■ , i^ 

Following the derivation given in Refs. [17|, [18[, where 
the non-relativistic limit is taken and the hadron wave 
functions are assumed to be those of a spherically sym- 
metric harmonie oscillator, the above equation can be 
reduced to 



^coal 



A N i TT S<Pvi<PkMki)f W (yi,ki) nnA 
9h 11 ~ 11 r m.. mT. sn.-\ ' ( 10 ) 



j=l "i í=\ 



J d 3 y i d 3 k i f l (k i 



where / (j/i, fcj) with yi and kí being, respectively, the 
internal (relative) spatial and momentům coordinates is 
the Wigner function associated with the internal (rela- 
tive) wave function. As in Ref. [13], we assume that the 
particles in the emission source are uniformly distributed 
in space and háve momentům distributions fj(j>j) given 
by the Boltzmann distribution of temperature T only for 
the transverse momentům Pj,Ti while the strong Bjorkcn 
correlation of cqual spatial (i]j) and momentům (Yj) ra- 
piditics is imposed for the longitudinal momentům, i.e., 



The coalescence model for particle production in nu- 
clear reactions is based on the sudden approximation by 
calculating the overlap of the density matrix of the con- 
stituents in an emission source with the Wigner func- 
tion of the produced particle [33||. The model has been 
extensively ušed to study light nuclei production in nu- 
clear reactions [34[ as well as hadron production from 
the quark-gluon plasma produced in relativistic heavy 
ion collisions |22h241 |35| . In the coalescence model, the 
number of hadrons of certain type h produced from the 



fj(Pj) ^(Yj -?7j)exp - 



2rrijT 



(11) 



with Tjj = log[(í 3 - + Zj)/(tj - Zj)]/2 and Yj = log[(Ej + 
Pj. z )/(Ej~pj, z )]/2 being the space-time and momentum- 
energy rapidities, respectively. From the relation 



7 ex P 



2rrijT 



exp 



P 2 

2MT 



\fiihl (12) 



with Pt and M denoting thc total transverse momentům 
and the total mass, respectively, we then obtain the 2- 
dimensional momentům distribution function of the con- 
stituents in Jacobi coordinates, fi(ki), 



fi(ki) <x exp 



A' 2 



2faT 



(13) 



where the reduced constituent masses fa are defined by 
1 1 1 



or explicitly 

Mi = 
M3 = 

l-H = 



fa m i+1 J2)=i 



?7iim2 m 3 (mi + m?) 

M2 ~ 



(14) 



mi + 7772 777l + "72 + "73 

7774(7771 +777 2 + 7Í7 3 ) 
777,1 + 7772 + 77^3 + 7774 ' 

7775(7771 + J772 + 777,3 + TfM) 



m 1 + Í77 2 + 7773 +7774+7775 



etc. (15) 



In the non-relativistic limit, the rapidity variables are 
simplified at midrapidities (Y = i] ~ 0) as r)j ~ Zj/tj and 
Yj ~ pj^ z /rrij. We can thus omit the contribution from 
thc longitudinal momentům in thc Wigner function / 
as long as the time where the coalescence takés pláce af- 
ter the collision is large compared with the internal time 
scale of the hadron, tj 3> 1/w, where u is the oscillator 
frequency. In this čase, the 3-dimensional momentům in- 
tegrations in Eq.(fT0|) reduces to 2-dimensional ones over 
fi(ki) and the Wigner functions in transverse momentům 
ki. The latter are given explicitly as 



faiVuh) =8 exp 



UJ 



i 2 2 



f^(m,h)^^8 + fafk^ 



W(m,ki) 



x exp 
16 " 



30 



VÍ 



.y±. _ 1.2 -.2 

2 K i "i 



20-^ + 15 - 20o?fc? + iafkf 



+I6yfkf - 8(77, • k t ) 



exp 



2 ^i ° i 7 ' 



(16) 



for the s-wave, p-wave, and cř-wave, respectively, with 
the parameters <7j = 1/^/fauJ related to the oscillator fre- 
quency w and the reduced constituent masses /i;. In Ap- 
pendix A, we derivc thc Wigner function for the d— wave 
statě from the harmonie oscillator wave functions. 

Carrying out the phase-space integrals in Eq. ([TO]) as 
shown in Appcndix B, we obtain the coalescence factor 



for each relative coordinate, 

_Jd 3 y t d 2 kJ l (k l )f w (y t ,k t 



F(ai,fa,k,T) 



fcPkiMkt) 
(4^ 2 ) 3 / 2 (2/i)!! 



l + 2faTa?{2h + iy. 



2faTaf 



IfaTaf 



(17) 



where ř, is the angular momentům of the wave function 
associated with the relative coordinate j/j. Combining 
these results, we obtain the following simple expression 
for the yield of hadrons from the coalescence model: 



jycoal 



N- 



N; 



3 = 1 3J 

(2ii)!l 



(47T<7?) 



2\3/2 



2faTaf 



(2ii + l)!] l(l + 2faTaf) 
gV{Mwf/ 2 (2T/uj) l 



(4tt) 3 / 2 



11 „.T//Vw.,W2 11 



(1 + 2T/uj) n + 

Tl-l 



?i+L-l 

(2/,)!! 



lí gjV (m J u J )^l\(2k + iyr 



(18) 



where žj is for an s-wave, 1 for a p-wave and 2 for a 
d-wave constituent, L = YH=i k, and M = X)"=i m í- 
Here we háve ušed the relation faof = 1/lu to convert 
the main dependence on li into the form of the orbital 
angular momentům sum L. For L > 2, the last factor 
in Eq. (fl"8]) depends on the way L is decomposed into k. 
For example, for L = 2 and n = 3, the combination 
(Zi,ř 2 ) = (1,1) gives a factor 4/9, while (h,l 2 ) = (2,0) 
leads to a factor 8/15. 



i. Quark coalescence 

To apply thc coalescence model to hadron production 
from the QGP at thc eritical temperature T c when the 
volume is Ve, we need to fix the oscillator frequency ap- 
propriately. This is doně by choosing the oscillator fre- 
quencies for light, strange, charmed, and bottom hadrons 
(w,o; s ,a; c and w&) in the quark coalescence to reproduce 
the yields of reference normál hadrons in the statistical 
model. These values are then ušed to predict the yields 
of exotic hadrons. 

For hadrons composed of light (up and down) quarks, 
we také the oscillator frequency u = 550 MeV to obtain 
in the coalescence model similar uj and p yields as in the 
statistical model as shown in Table. IIIII 

For hadrons composed of light and strange quarks, the 
parameter w s in thc coalescence model is determined by 
fitting the statistical model prediction for A(1115) in- 
cluding thc contribution from resonance decays. Taking 



TABLE III: Yields of normál hadrons at RHIC and LHC 
in the coalescence and statistical models with oscillator řre- 
quencies co = 550 MeV,cj s = 519 MeV,o; c = 385 MeV, and 
uib = 338 MeV are determined by fitting the statistical model 
results for A(1115), A c (2286), and Af,(5620) marked with (*) 
at RHIC after taking account of resonance decays. Numbers 
in the parentheses are those without the decay contribution. 





particle 


RHIC 


LHC 




coal. stát. 


coal. 


stát. 




w(782) 


44.2 40.2 


119 


108 


qq 


p(770) 
K* (892) 


132 127 
41.2 47.2 


358 
111 


342 
135 




A"* (892) 


41.2 52.9 


111 


135 


qqs 


A(1115) 
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into account states in the octet and decuplet reprcscnta- 
tions that decay dominantly to A(1115), we obtain the 
following result for heavy ion collisions at RHIC: 



N 
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A(1115) 



yvrstat 

iV A(1115) 



Arstat 
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/vrstat 
iV ~(1318) 
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i Arstat 
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Arstat 
iv n-(1672) 

x 13.57 + 4.73 



0.87- 



0.11 



x 10.91 + 3.42 



29. 



0.81 

(19) 



We will use this parameter to estimate the yield of other 
hadrons that are composed of light quarks and strange 
quarks. As shown in Table ITTl| this value leads to a yield 
of A(1520), which has the s quark in the p-wave statě, 
in the coalescence model that is smaller than that in the 
statistical model as hrst pointed out in Ref. |31| . 

The oscillator frequency for charmed hadrons is fixcd 
to reproduce the A c (2286) yield including the fecd-down 
contribution [36j but without taking into considcration 
the effect of the diquarks [Š7J]. For the A c (2286) yield, we 
consider only the contribution from £ c (2455), £ c (2520) 
and A c (2625) decays as states of higher masses are neg- 
ligible, that is, 
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A c (2286) 
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Arstat i /vrstat 

JV S C (2455) + JV E C (2520) 



+0.67 X ^2625) 



0.139 
0.602 



0.177 + 0.254 



0.67 x 0.048 

(21) 



at RHIC. Fitting again the total yield of A c (2286) calcu- 
lated in the statistical model to that in the coalescence 
model for the samé resonances as given in Eq. (|21[) , 



=0.602, 



0.173 + 0.346 + 0.67x0.037 



(22) 



we obtain tu c = 385 MeV for the charm quark mass m c = 
1500 MeV. 

The oscillator frequency for bottom hadrons u>b = 
338 MeV is obtained by fitting the sum of the statisti- 
cal model results for A;, (5620) and the contribution from 
£ b (5810) and £*(5830) decays at RHIC, 
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-3 



1.73 x 10" 3 

(23) 
3 + 2.170 x 10~ 3 

(24) 



In the above formula, 0.87 and 0.11/3 in the bracket rep- 
resent, respectively, the branching ratios of £(1385) — > 
A + 7T and £(1385) ->• £° + tt in the £(1385) decay. All 
numbers are calculated at Th and Vh with /i s — 10 AleV 
and hb = 20 MeV for RHIC. To reproduce the total yield 
within the coalescence model with the constituent quark 
masses m u .d = 300 MeV and m s = 500 MeV, we need 
uj s = 519 MeV after taking into account the samé feed- 
down contributions as in Eq. (|19j). Spccifically, we háve 
from the coalescence model 



^Aaiísf 1 =3 - 01 + ~ * 



3 

0.87- 



0.11 



=29.8. 



2.20 

x 18.07 + 4.40 + 0.78 

(20) 



with the bottom quark mass m& = 4700 MeV. 

Since the oscillator frequencies are related to the sizes 
of hadrons [13, 13 ; samé values as determined at RHIC 
are ušed in the coalescence calculations for heavy ion col- 
lisions at LHC. Using the samé ui values as those for nor- 
mál hadrons, we then see from Eq. (fT8|) that the addition 
of a s-wave, p-wave, or d-wave quark leads to, respec- 
tively, a coalescence factor 



1 Nj (4^t| ) 3 / 2 
J l V(í + 2n i Taf) 
1 Ni 2 (47ra i 2 ) 3 / 2 2/z z T ( 7f 
9l V 3 (1 + 2/ííTct?) 2 
1 N 8 (4 7 r ť T 2 ) 3 / 2 (2^Ta 2 ) 2 
9l V 15 (1 + 2 AÍJ Ta 2 ) 3 



-0.360 



-0.093 



-0.029. 



(25) 



The production of multiquark hadrons involves more s- 
, p- and d-wave coalescence factors and is hence gen- 
erally suppressed. Moreover, the d-wave coalescence is 
supprcsscd in comparison with the p-wave coalescence, 
which is further suppressed relative to the s— wave coa- 
lescence 13111 . 



2. Hadron coalescence 

For the yields of wcakly bound hadronic molecules 
from the coalescence of hadrons, they are evaluated at 
the kinetic freeze-out temperature Tp and volume Vf- 
The oscillator frcquencies needed for hadronic molecules 
in the hadron coalescence is related to their mean square 
distances (r 2 ) between the two constituent hadrons. For 
a hadronic molccule in the relative s-wave statě, the os- 
cillator frequency is given by 



2/ii?(r 2 



(26) 



where /j,r = m\mil(rn\ + m-i) is the reduce mass. The 
mean square distance of the hadronic molecule can be 
further related to its binding energy B via the scattering 
length a,Q of the two interacting constituent hadrons, i.e., 



^oí2 for the diomega (ížíž) + predicted by the chiral 
quark model [38| can be calculated from Eq. (J2GI) . 



Wfif2 — 



2//nn('" 2 )nn 
3 197.3 2 



2 1672.45/2 x 0.84 2 



98.8 MeV (29) 



where \Z(r 2 ) nn = 0.84 fm [38(] has been ušed. By the 
samé token, wc can calculate the oscillator frequencies 
for all the hadronic molecules, and the results are sum- 
marized in Table IPV1 



III. EXOTIC HADRONS 

In this Section, we briefly discuss the properties, such 
as the quantum numbers and possible decay modes, and 
the current theoretical and experimental status of the 
exotic hadrons includcd in the present study. For conve- 
nience, we classify these hadrons into exotic mesons (Sub- 
scc. A), exotic baryons (Subsec. B) and exotic dibaryons 
(Subsec. C) as shown in Table IIV1 

A. Exotic mesons 



D 



'2fj, R al 



« = a , 



(27) 



These relations are valid when the binding energy is small 
and the scattering length is large compared to the range 
of the hadronic interaction, and they can be easily ob- 
tained as follows. Using the relation kcotSo(k) = —1/ao 
for the relative momentům k — > between the s-wave 
scattering phase shift Sq and the scattering length ao, 
the S matrix for two interacting hadrons at low energies 
can be approximated as 



O 2io (fc) 



-i + ife 

ao 

— ik 



j_ 

ao 



(28) 



which has a pole at k = i/a 0l corresponding to a bound 
statě with the binding energy given by the hrst equation 
in Eq. (|27|) and the radiál wave function outside the inter- 
action range Ub(r) ~ e~ r > ao . Assuming that ao is much 
larger than the interaction range and using the above 
form of the wave function for the whole region, we obtain 
the mean square distance given by the second equation 
in Eq. (|2"?| . For a wcakly bound two-body states, we 
thus obtain from Eqs. (f2"6")l and (J2"T)) the simplc relation 
uj = QB. We notě that (r 2 ) is the mean square distance 
in the relative coordinate, and it is not the squared mean 
rádius from the center-of-mass. 

For example, the oscillator frequency for /o (980) can 



be obtained from il^ (9 80 ) = 6 x B 



-5/0(980) — M* 



■ +Mft-M 



/o (980) 



67.8 MeV with 



/o (980) 



493.7+497.6-980: 



11.3 MeV. As another example, the oscillator frequency 



For exotic mesons, we include the following: 

1. /o (980): This 1 = scalar particle together with 
the 1=1 eto(980) are members of the scalar nonet 
that has been thought to be composed of multi- 
quark configurations [0, uj ■ If /o (980) is a mem- 
ber of the multiquark configurations, its wave func- 
tion then has a hiddcn šs component. Assuming 
that they are composed of quark-antiquark pair, 
their wave functions would be /o (980) ~ sš and 
ao(980) ~ (uu — dď)/y/2. An early QCD sum rule 
analysis suggested, on the other hand, that /o (980) 
and ao(980) were just the 1 = and 1 combinations 
of uu and dď [39| . There are also modcls in which 
/o (980) is a KK molecule. While there seems to be 
consensus from lattice calculations that this parti- 
cle is a tetraquark statě, the situation is not at all 
clear because of the treatment of disconnected dia- 
grams [40l |. 

2. if (1460): This is an exited statě of kaon with 
J F = - on the particle list of the Particle Data 
Group (PDG) but is omitted from the summary 
table [41| . It was observed in the partial wave 
analysis of the Kirir finál statě in elementary re- 
actions. Recently, this resonance has been stud- 
ied theoretically in the context of the meson dy- 
namics. In Ref. [43] this kaon was obtained from 
the K- /o (980) s-wave two-body dynamics with the 
/o (980) dynamically generated in the KK and irir 
coupled channels. A later study in Ref. [43( con- 
sidered the three-body couplcd-channel dynamics 



TABLE IV: List of exotic hadrons discussed in this páper. Shown are the mass (m), degeneracy (g), isospin (I), spin and parity 
(J p ), the quark structure (2q/3/q/6q and Aq/5q/Sq), molecular configuration (Mol.) and corresponding oscillator frequency 
(wmol), and decay mode of a hadron. For the WMoi., it is fixed by the binding energies B of hadrons (oj ~ 6 x B, marked (B)) 
or their mean square distances (r 2 ) (u) ~ 3/2/x (r 2 ), marked (R)). In the čase of three-body molecular configurations for exotic 
dibaryons, we adopt the WMoi. as that for the subsystém, as marked (T). Further marked by *' are undetermined quantum 
numbers of existing particles, by " particles which are not yet established, and by ' particles which are newly predicted by 
theoretical models. 
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of KKK in a Faddeev approach and found a very 
similar kaonic excitation. A non-relativistic poten- 
tial model for the KKK systém was also ušed in 
Ref. [43| , and a quasi-bound statě of binding energy 
21 MeV and root mean square rádius 1.6 fm was 
obtained. Thus, the AT(1460) could be understood 
as a KKK hadronic molecular statě. 

3. D s j(2317): This statě was first observed by the 
BaBar Collaboration [2| through the Dftr channel 
in inclusivc e+e - annihilation. Its measured mass 
is approximatcly 160 MeV bclow the prediction of 
the verv successful quark model for the charmed 
meson [44|. Due to its low mass, the structure of 
D s j(2317) has been under extensive debatě. It has 
been interpreted as a cš statě [4a - [49l ]. two-meson 
molecular statě I50l . [5I| , D — K mixing [52| , four- 
quark statě [53I 456J qr a mixture of two-meson and 
four-quark states [57J | . 



T^ c (udčč): The set of tetraquarks with two heavy 
quarks wcre first considcred in Ref. [58( . The struc- 
ture with [ ud] diquark is expected to be particu- 
larly stable [53, |60] and could be bound against the 
strong decay into D\D. The quantum number is 
J p = 1+ with 1 = 0; hence the decay into DD is 
forbiddcn due to the angular momentům conserva- 
tion. Estimates based on the simple color-spin in- 
teraction suggest the mass to be 3796 MeV |13l[l4|. 
The hadronic decay mode of T cc is D*~D° if its 
mass is above the threshold and Ď°D°it~ if be- 
low. If the T cc is strongly bound, it can then decay 
weakly to D*~K + tt~ with a lifetime similar to that 
of the Ď meson. A molecular statě with the samé 
quantum number was also predicted to exist within 



the pion-exchangc model |61| . The production of 
doubly charmed hadrons in heavy ion collisions at 
RHIC was discussed in Ref. K32| , and also estimated 



for LHC [13( . This statě could also be searched for 
at Belle. In particular, the search should be simi- 
lar to that for the doubly charmed baryon such as 
the S cc . Here are the two possibilities. The first 
one is through the decay of the B meson. Unfor- 
tunately the search at Belle was not successful (63| 
as the dominant subprocess was the weak decay 
of the b quark into a č by cmitting a W + (— > cš). 
Therefore, while the cč is produced, the cc creation 
might be highly suppressed. On the other hand, 
another more feasible search is in the continuum 
background where two cč pairs are known to be 
produced in the reaction e + e~ — > J/^X(3940) [64j | . 

A (3872): The Belle collaboration found this par- 
ticle in the B+ -> X(3872)K+ -> J / ip-K+ -k~ K + 
decay Q. CDF, DO, and BaBar háve confirmed its 
existence, and the current world average mass is 
3871.2 ±0.39 MeV. Although the new BaBar result 
favors the J PC = 2 h assignment |65|, the estab- 
lishcd properties of the X(3872) are in confhct with 
this assignment [66|, HŤj]. Therefore, the favored 
quantum numbers are J PC = 1 ++ with the isospin 
violating decay modes. This particle was predicted 
in Ref. [68[ as a DĎ*° bound statě within a meson- 
exchange model. It was shown in this study that in 
the DĎ* sector the one-pion exchange interaction 
alone is strong enough to form a molecular statě 
that is bound by approximately 50 MeV. In this 
čase, other molecular states of D*Ď* , D\D* , DiĎ 
and DqD are also expected to exist as the pion ex- 
change is allowed in these channels as well. In other 
studies, it has been claimed that the X(3872) has 
the admixture of the čc component and is thus a 
tetraquark hadron [5j. Its dominant decay modes 



include J/ijjtt' 



J/^tt+tt", andD a Dir 



6. Z + (4:4:30): The Belle collaboration observed this 
charged statě in B + — s- KiP'tt + through its decay 
into ý/TT + Q. The reported mass and width are 
M = 4433 MeV and T = 45±^|±fg MeV, respec- 
tivcly. The reported mass is closc to the D\Ď* 
threshold and hence the possible structure for this 
statě is cither a molecular or a tetraquark statě 
Q- As commented in the discussion on A (3872), 
the one-pion exchange interaction could bound a 
D\D* molecular statě. So far, no other experiment 
has confirmed this finding. In particular, BaBar 
[69| also searched the Z" (4430) signatuře in four 
decay modes: B — > iJm~K, where ip = J/ip or 
ip' and K — Kg or K + . No significant evidence 
for a signál peak was found in any of the investi- 
gated processes. After the failure of the BaBar col- 
laboration in confirming the Z~(4430) mass peak, 
Belle has performed a reanalysis of their data that 
took detailed account of possible reflections from 
the Kir~ channcl. From a full Dalitz plot reanaly- 
sis of their data, Belle has confirmed the observa- 
tion of the Z + (4430) signál with a 6.4er peak sig- 



nificancc. The updated Z + (4430) parameters are: 
M = (4433t^+i9)MeV and T = (109Ígí^)MeV 
[Ť0|. If confirmed, the Z+(4430) is the first prime 
candidate for an exotic particle. Considering the 
Z + (4430) as a loosely bound s-wave D\D* molecu- 
lar statě, the allowed angular momentům and par- 
ity are J p = 0~,1~,2~, although the 2~ assign- 
ment is probably suppressed in the B + — > Z + K 
decay by the small phase spacc. Among the re- 
maining possible _ and 1~ states, the former will 
be more stable as the latter can also decay to D\Ď 
in s-wave. Hence the 0~ quantum number is fa- 
vored. 

Věry recently the BELLE collaboration [7l| re- 
ported the observation of two narrow charged 
structures in the hidden-bottom decay channels 
■n ±r l{nS) (n = 1,2,3) and n ± h b (mP) (m = 1,2) 
of T(55). The measured masses (widths) of these 
two structures are, in units of MeV, Mz b = 10610 
(T Zb = 15.6 ± 2.5) and M z > b = 10650 (F z - = 
14.4 ± 3.2), respectively. The analysis of the Z\, 
states decay in the channel Z^ — > T(2S')7r ± fa- 
vors the I (J ) = 1 + (1 + ) assignment. Sincc 
the masses of these two states are very close to 
the BB* (10604.6 MeV) and B*B* (10650.2 MeV) 
thresholds, they are ideál candidates for these 
molecular states. It is interesting to notice that the 
decay channels of the Z b states are very similar to 
the decay channel of Z + (4430). 

7. T® b (udčb): Both the light diquark and heavy anti- 
diquark are scalar diquarks so that J p = + with 
1 = 0. This particle could be strongly bound with 
a mass of 7149 MeV[H|. As it has the Ď°B° com- 
ponent, it can decay weakly via T® b — > K + it~ + 
K+tv-. 



B. Exotic baryons 

For exotic baryons, we consider the following: 

1. A(1405): This resonant statě with / = 0, J p = 
1/2- , mass 1406±4 MeV, and width 50±2 MeV [U 
has been considered as a quasi-bound statě of the 
KN systém [9j, even before the establishment of 
the QCD. The modern theoretical approach based 
on the chiral dynamics within the unit ary frame- 
work (the chiral unitary approach) |l(| [HI, I72h76J 
also suggests that this resonant statě is dynamically 
generated in the meson-baryon scattering includ- 
ing the KN and 7rS channels, and is dominated by 
the meson-baryon molecular component [77] . The 
mean square distance between K and N in the 
A(1405) is evaluated to be (r 2 ) = 2.7 fm 2 in the 
chiral unitary approach [Ť8|], in which the A(1405) 
peak appears around 1420 MeV in the KN channel. 
This value leads to a smallcr oscillator frequeney 



10 



for the bound statě (wmoI. = 20.5 MeV) than that 
fixed by the binding energy with the A(1405) mass 
of 1405 MeV. Its dominant decay mode is 7r£ in the 
1 = channcl, but there may be the possibility of 
observing it in the 7A decay mode. 

2. <d + (uuddš): This fiavor exotic baryon with 
strangeness S = +1, J p = 1/2+ and 1 = was 
predicted in the chiral soliton model [79|. The in- 
triguing features are its light mass of 1540 MeV 
and narrow width, which partly motivated the first 
cxperimental observation by LEPS [12| ■ Although 
the LEPS result was confirmed by sever al other 
cxperimcnts, it was subscquently followed by neg- 
ative results with high statistics such as in the 
high energy collision experiment by PHENIX (801 
and the low energy photoproduction experiment by 
CLAS [H|]. On the other hand, the recent LEPS 
result maintains a positive signál [82| . The most 
suitablc hadronic decay mode for its identification 
in an inclusive experiment is 0+ — > K°p. The pro- 
duction rate in heavy ion collisions has been esti- 
mated in the statistical model |83l - l85| and in the 
coalescence model [131. ^ n the quark model, the 
spin-parity of 0+ is 1/2" if all five quarks are in the 
s-wave orbit but is 1/2+ after including the strong 
diquark correlation (86[. The possibility of a 3/2~ 
statě has also been proposed to explain its narrow 
width [87J. A recent comprehensive QCD sum rule 
study favors, however, the 3/2+ assignment |88j . 

3. KKN: The quasibound statě of KKN was pre- 
dicted to be a hadronic molecular statě of N* with 
a mass of 1910 MeV, J p = 1/2+ and / = 1/2 in 
the variational calculation using a hadronic two- 
body potential (Šíj. This statě was confirmed by 
a coupled-channels Faddeev calculation [90J . It 
has been interpreted as a coexistence statě of the 
A(1405)AT and ao(980)V clusters, and its main de- 
cay modes are thus the KirH from the A(1405) de- 
cay and the irr/N from the ao (980) decay. Since the 
KKN is a hadronic molecular statě, it has a large 
spatial distribution. The root mean square rádius 
is found to be 1.7 fm and the interhadron distances 
are larger than 2 fm [89j |. 

4. ĎN(Ď*N) and BN(B*N): The quark contents of 
these hadron molecular states are similar to the 
c (b) but with the different quantum numbers of 
J = 1/2" and 1 = 0. Recently, based on the 
heavy quark symmetry, a pion induced bound DN- 
Ď*N molecular statě was predicted to exist with a 
binding en ergy of a few MeV below the threshold of 
2806 MeV [T5[ . This is a shallow bound statě com- 
pared to the deeply bound 9 C of about 100 MeV 
below the threshold. An easily identifiable decay 
mode is K + iv~ / K~p. The BN molecule would be 
more stable, because the heavy quark symmetry 
amplifies the strong mixing bctwcen BN and B* N 



and thus suppresses the kinetic energy. The mass 
of the BN molecule was predicted to be a few tens 
MeV below the threshold of 6218 MeV, and the 
possible weak decay mode is K + tt~tt~ + 7r+ + p. 
In a more recent study [ljl, the ĎN (BN) was 
also predicted to háve a resonance statě between 
the ĎN (BN) and Ď*N (B*N) thresholds, with 
the mass 2929 (6226) MeV and quantum numbers 
J p = 3/2~ and 1 = 0. Similar to the ĎN (BN) 
bound states, this resonance statě is induced by 
the pion exchange within the heavy quark symme- 
try. This resonance statě can also be regarded as 
a bound statě of Ď* (B*) and N with respect to 
the Ď*N (B*N) threshold like that of the Fesh- 
bach resonance. Its decay width of 19 (0.12) MeV 
to the ĎN (BN) via the strong interaction is very 
narrow as a result of the suppression by the d-wave 
centrifugal barrier in the finál ĎN (BN) statě. We 
háve fixed the oscillator frequeney lo from the root 
mean square distance, \/ (r 2 ) = 3.8(1.7) fm for the 
ĎN (BN) molecular statě [Tí|, and the samé u> is 
ušed for its excited statě, Ď*N (B*N). 

5. @ c (uuddč) and Qb(uuddb): The bound Skyrmion 
approach predicted the bound exotic hadron O c 
with the mass of 2650 MeV and quantum num- 
bers J p = 1/2+ and 7 = J9l|. There was one 
experiment reporting a positive signál [92| at a 
mass around 3.1 GeV, but 110 confirmation exists 
so far |23|. An easily identifiable decay mode is 
K + Tr~ir~p if the statě is strongly bound and D*~p 
if it is a resonant statě. Similarly, the Ob mass was 
predicted to be 5207 MeV with the samé quantum 
numbers J p = 1/2+ and 1 = 0. The possible weak 
decay mode is K + Tr~n~ + 7r+ + p. 

6. O cs (udusč): In the quark model including the 
color-spin interaction, the J p = 1/2 - and 1 = 
five-quark statě can be bound and becomes sta- 
ble against the strong decay [9J, |95|. The mass 
is predicted to be 2920-2930 MeV, depending on 
the model parameters. This statě was searched for 
in the Fermilab E791 experiment through the ({mp 
mode [9a] and the K*°K~p mode [97f, and the re- 
sults are so far negative. 



C. Exotic dibaryons 

The exotic dibaryons included in the present study are: 

1. H dibaryon: This particle was first predicted in 
Rcf. [8j as a deeply bound statě below the AA 
threshold. Dcspitc extensive searches such as in 
the BNL-E885 experiments [98J], deeply bound H 
dibaryons háve not been observed. The discovery 
of the double-A hypernucleus ( AA He) in the Nagara 
event [99| finally exeluded the possibility of deeply 
bound H dibaryon, as the two A particles can decay 



11 



strongly into the core nucleus and the H dibaryon if 
the H mass is below the AA threshold by more than 
the AA separation energy (B AA = 7.25 ± 0.19Ío'.u 
MeV) j99[ . As a result, we now háve only a 
narrow window for the H particle to be bound, 
< Bh ;$ 7 MeV, where Bh is the binding en- 
ergy of H from the AA threshold. The reason why 
the H particle docs not strongly bound may be due 
to the instanton-induced determinant (Kobayashi- 
Maskawa-'t Hooft) |l00j interaction, which acts re- 
pulsively in the H channel and may cancel the 
strong color-spin attraction, as demonstrated in the 
quark-cluster model |lOl| . There is, however, still a 
possibility that the H particle exists as a resonance. 
In the KEK-E522 experiment, an enhancement in 
the A A invariant mass spectr um is observed at 10- 
20 MeV above the threshold [lO^, while the signif- 
icance as a peak is only around 2a. This peak-like 
enhanceme nt can not be explained by the final-state 
interaction |l03J ] . Recent lattice calculations háve 
suggested th at a bound statě pole exists around the 
SU(3) limit [l0|- With the realistic SU(3) break- 
ing, this bound statě pole would be shifted to a 
weakly bound statě or a resonant statě between the 
A A and SiV thresholds. If we apply the low energy 
scattering formula (Eq. (j2"7| ). the rms rádius of the 
H resonance as a bound statě of EiV [l02| may be 
evaluated to be 0.9 — 1.3 fm. 

2. KNN: Motivated by the existence of the A(1405) 
resonance below the KN threshold, the possibil- 
ity o f bou nd _ftT-nuclear systems was proposed in 
Ref. J105J | based on a phenomenological KN po- 
tential. Since then the simplest KNN systém 
has been intensively studicd both thcorctically and 
cxpcrimcntally. While the experiment by FIN- 
UDA [l06J indicates a peak structure at 2255 MeV 
in the KN invariant mass spectrum, the interpre- 
tation of the p eak as the KNN statě is still con- 
troversial jl07| . Recent rigorous few-body calcula- 
tions for the KNN systém indicate that the sys- 
tém bounds in the J p — _ and 7 = 1/2 chan- 
nel [l08l - llll | . With a suitable treatment of the en- 
ergy dependence of the KN interaction, the mass 
of the KNN sy stém is predicted to be about 2350 
MeV |l!2l4ll4 |. In heavy ion collisions, this statě 
can be observed in the KN or irllN invariant mass 
spectrum. 

3. (ížíž)o+: This is a deeply bound six- quark statě pre- 
dicted by the chiral quark model |38|, Ill5| . It has a 
large binding energy of about 116 MeV and a small 
root mean square distance of 0.84 fm between the 
two fžs. Because of its large strangeness content, 
it is stable against strong hadronic decays and pos- 
sesses the weak decays (fžfž) + — > 7r~ +5° + fž~ and 
(f2f2)o+ — > Ti" + S~ + Í2~ with a mean lifetime esti- 
mated to be about four times longer than the free 
íl lifetime of 0.822 x 10~ 10 sec. Apart from these 



conventional decay modes, the nonmesonic decay 
(OO) + — > 5~ + íž _ is also possible; and the esti- 
mated lifetime of (íifi) + for this process is twice 
the free íž lifetime. Thus, instead of direct observa- 
tion, the (ížíž) + may also be detected in the S~0~ 
invariant mass distribution. 

4. H++: This dibarvon with J p = 0° and I = 1 is 
predicted in Ref. J14j . It is expected to be strongly 
bound as it is composed of [ud] , [us] , and [uc] scalar 
diquarks, and one of which has to break in order for 
H+ + to fall apart and to decay to p + S + . The life- 
time is expected to be similar to that of S + , and the 
dominant hadronic weak decay mode is expected to 
bc p + AK~tv + tt + . 

5. ĎNN and BNN: The attractive ĎN and BN in- 
teractions would lead to bound D and B mesons 
in the nuclear medium as their binding energies in- 
crease with inereasing nuclear density. The DNN 
and BNN molecular states predicted in Ref. [15( 
are thus nuclei with minimum baryon number. The 
quantum numbers can be J p = 0~ or 1 _ and 
1 = with different types of weak decay mode. The 
ĎNN states háve the decay modes K + ir~ir~ +p+p 
for _ and K + tt~ + d for 1~ with all charged par- 
ticles in the hnal statě, while the BNN states háve 
the decay modes K + tt~ + tt + + p + p for 0~ and 
K + ir~ + n + + d for 1~. Thercfore, the experimen- 
tal observation of these decays makes it possible 
to determine both the spins and parities of these 
hadronic molecular states. 



IV. YIELDS OF EXOTIC HADRONS IN HEAVY 
ION COLLISIONS 

We show in this Section the expected yields of exotic 
hadrons deseribed in the previous Section from centrál 
Au+Au collisions at v / š]v r žv = 200 GeV at RHIC and cen- 
trál Pb+Pb collisions at ^š^T = 5.5 TeV at LHC. They 
includc rcsults for all possible structure configurations, 
e.g. multiquark hadrons and hadronic molecules, calcu- 
latcd from the coalescence model in addition to those 
estimated from the statistical model. These results are 
shown in Table [V] We also give some discussions on the 
obtained rcsults. 

Comparisons of the yields in the 2q/3q/6q column to 
those in the 4q/5q/8q column in Table [V] show that for 
most of the hadronic states considered here, the yield 
from the coalescence model for the compact multiquark 
statě is smaller than that for the usual quark configura- 
tion as a result of the suppression duc to the coalescence 
of additional quarks indicated in Eq. (f2~5"j) . For the samé 
statě, the yield from the coalescence model for the molec- 
ular configuration is, however, larger than that from the 
statistical model prediction as seen from comparing the 
yields in the Mol. column to those in the Stát. column 
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TABLE V: Exotic hadron yields in centrál Au+Au collisions at y^šjviv = 200 GeV at RHIC and in centrál Pb+Pb collisions at 
,/snn — 5.5 TeV at LHC from the quark coalescence (2q/3q/6q and 4q/5q/8q) and the hadron coalescence (Mol.) as well as 
from the statistical model (Stát.) 





RHIC 


LHC 




2^/39/6^ 


4q/5q/8q 


Mol. 


Stát. 


2q/3q/6q 


4g/5g/8g 


Mol. 


Stát. 


Mesons 
/o (980) 
a (980) 
K{U60) 
D s (2317) 
T 1 t) 

-*■ CC 

X(3872) 

Z + (4430) t) 

T o t) 

M cb 


3.8, 0.73(sš) 
11 

1.3 x 10" 2 

1.0 x 10"* 


0.10 

0.31 

0.59 
2.1 x 10" 3 
4.0 x 10" 5 

4.0 x 10" 5 
1.3 x 10" 5 

6.1 x 10" 8 


13 

40 

3.6 
1.6 x 10" 2 
2.4 x 10" 5 
7.8 x 10"* 
2.0 x 10" 5 
1.8 x 10" 7 


5.6 

17 

1.3 

5.6 x 10" 2 

4.3 x 10~* 
2.9 x 10~* 

1.4 x 10" 5 
6.9 x 10" 7 


10, 2.0 (ss) 
31 

8.7 x 10" 2 

1.7 x 10" 3 


0.28 

0.83 

1.6 
1.4 x 10" 2 
6.6 x 10"* 
6.6 x 10"* 
2.1 x 10"* 
6.1 x 10" 6 


36 

1.1 x 10 2 

9.3 

0.10 

4.1 x 10"* 

1.3 x 10" 2 

3.4 x 10"* 
1.9 x 10" 5 


15 

46 

3.2 

0.35 

7.1 x 10" 3 

4.7 x 10" 3 
2.4 x 10"* 

6.8 x 10" 5 


Baryons 
A(1405) 

e +t) 

KKN ť 
DN t) 

Ď*N f) 

BN f) 
B*N t) 


0.81 


0.11 
2.9 x 10" 2 
1.9 x 10" 2 
2.9 x 10" 3 
7.1 x 10"* 
5.9 x 10"* 
1.9 x 10" 5 
5.3 x 10" 6 


1.8-8.3 

1.7 
4.6 x 10" 2 
4.5 x 10" 2 

8.0 x 10" 5 
1.2 x 10"* 


1.7 
1.0 

0.28 
1.0 x 10" 2 
1.0 x 10" 2 
7.2 x 10" 3 
3.9 x 10" 5 
6.6 x 10" 5 


2.2 


0.29 

7.8 x 10" 2 
5.2 x 10" 2 

2.0 x 10" 2 
4.7 x 10" 3 

3.9 x 10" 3 
7.7 x 10"* 

2.1 x 10"* 


4.7-21 

4.2 
0.28 
0.27 

2.8 x 10" 3 
4.4 x 10" 3 


4.2 
2.3 
0.67 

6.1 x 10" 2 

6.2 x 10" 2 
4.5 x 10" 2 
1.4 x 10" 3 
2.4 x 10" 3 


Dibaryons 

Jí" 

if JVAr f) 

fifi t) 

//++ t) 

dnnv 

BNN t) 


3.0 x 10" 3 
5.0 x 10" 3 
3.2 x 10" 5 
3.0 x 10"* 


5.1 x 10"* 

2.9 x 10" 5 
2.3 x 10" 7 


1.6 x 10" 2 
0.011-0.24 
1.5 x 10" 5 
3.3 x 10"* 
1.8 x 10" 3 
1.2 x 10" 6 


1.3 x 10" 2 
1.6 x 10" 2 

6.4 x 10" 5 

7.5 x 10~* 
7.9 x 10" 5 
2.4 x 10" 7 


8.2 x 10" 3 

1.3 x 10" 2 
8.6 x 10" 5 
2.0 x 10" 3 


1.4 x 10" 3 

2.0 x 10"* 
9.2 x 10" 6 


3.8 x 10" 2 
0.026 - 0.54 

4.4 x 10" 5 

1.9 x 10" 3 
9.8 x 10" 3 
3.7 x 10" 5 


3.2 x 10" 2 
3.7 x 10" 2 
1.9 x 10"* 
4.2 x 10~ 3 
4.2 x 10"* 
7.6 x 10" 6 



in Table [V] This is in contrast to high energy pp colli- 
sions, where molecular configurations with s mall b inding 
cnergies are hard to be produced at high pt |116| . 

To see more clcarly the cffcct of the structure of a 
hadron on its production in heavy ion collisions, wc show 
in Fig. [1] the ratio of the coalescence model results to 
those from the statistical model, 



Rť = Nr l /N s h 



(30) 



Generally, the ratios for the Iq and 3q configurations are 
within the range of 0.2 < Rh < 2 for normál hadrons, 
which is shown by the grcy zonc. This also applies to 
exotic hadrons with most of the 2q/Zq configurations, 
as shown by the triangles. We observe that the coales- 
cence yield for an exotic multiquark hadron (diamond) is 
smaller than those for the usual quark configuration and 
from the statistical model predictions. This is consistent 
with the naivě expectation that the multiquark coales- 
cence becomes suppressed as the quark number inereases. 



The tetraquark states /o (980) and ao(980) are typical ex- 
amples. In these hadrons, tetraquark configurations nec- 
cssarily involve strange quarks, and they are thus more 
suppressed. This suppression also applies to the 5g states 
in exotic baryons (A(1405), e+(1530), KKN, and © cs ) 
and the 8q statě in the KNN. 

The yields of hadronic moleculars from the hadron coa- 
lescence (circles) depend strongly on the sizes of hadrons. 
For deeply bound and compact hadronic molecules, their 
yields are comparablc to or smaller than the predictions 
of the statistical model. On the other hand, the loosely 
bound extended molecules would be formed abundantly. 
One typical example is A(1405) that are shown in Ta- 
ble HV] for the two cases of u = 20.5 and 174 MeV. The 
smaller value, corresponding to a small binding |75l.[76{. 
gives a larger size for the A(1405). As a result, antikaons 
produced in heavy ion collisions háve larger probabili- 
ties of coalescing with nucleons, resulting in more abun- 
dant production of A(1405). On the other hand, for the 
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Coalescence / Statistical model ratio at LHC 
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FIG. 1: (Color online) Ratio of the yield of an exotic hadron in the coalescence model to that of the statistical model. 



larger u value, which corresponds to the čase that the 
pole position of the S-matrix for the two-hadron interac- 
tion is around 1405 MeV, the A(1405) can be regardcd 
as a deeply bound statě and thus has a smaller size, and 
hence its yield becomes smaller. 



V. DISCUSSIONS 

Our results based on the coalescence model for hadron 
production in relativistic heavy ion collisions háve indi- 
cated that their yields are strongly dependent on their 
structures. Therefore, measuring the yields of exotic 
hadrons allows us to infer the internal configuration of 
exotic hadrons [lj| [20( ■ For example, we háve mentioned 
in Sec. IIVI A that as possible configurations for /o (980), 
quark-antiquark pairs (~ sš, uu, and dd), a tetraquark 
statě, and a KK hadronic molecule háve been proposed. 
To conhrm its structure, wc refer to preliminary data 
from the STAR Collaboration for the production yield 
ratios of /o (980), 7r, and p° |ll7j . From these results we 
find that the measured yield of /o(980) is close to 8, which 
means that it is more probable for the /o (980) to be pro- 
duced as a hadronic molecule statě than a tetraquark 
statě (See the order of magnitudc difference between the 
yield in the Aq/5q/8q column and that in the Mol. col- 



umn in Table fVT) . Therefore, we conclude that the STAR 
data seem to rule out a dominant tetraquark configu- 
ration for the /o (980). Further experimental efforts to 
reduce the error bar are thus highly dcsirablc. 

For some exotic hadrons, our results show that the 
yields are similar for the hadronic and the molecular 
configuration, despite the difference in the coalescence 
temperatures Tq and Tp. This can be attributed to 
the larger size of the molecular configuration. Assum- 
ing other factors are similar, the s-wave factors involved 
in the coalescence at Tp are similar to those at Tc as long 
as the relevant molecular size is related to the hadron size 
as (7(7 = {Vc/VpY^ap as can be inferred from Eq. ([25]) 
after neglecting the temperature dependence in the de- 
nominator. If wc additionally assume that the volume 
scales as V cx T -3 , we find that the condition for the 
molecular coalescence to be similar to two-quark coales- 
cence is that the molecular size scales as <jp = acTc/Tp, 
which is more or less satisfied by some exotic hadrons 
considered here, such as the D s (2317). 

Our study also shows the interesting result that the 
ratio R^ s of the yield in the coalescence model to that 
in the statistical model is almost the samé at RHIC and 
LHC. This similarity comes from the univcrsal fcature of 
the QCD phasc transition; the common eritical temper- 
ature and the common volume ratio Vc/Vh- In the non- 
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FIG. 2: (Color online) Ratio of hadron yields at RHIC in the coalescence model to those in the statistical model as a function 
of mass. 



relativistic approximation shown in Eq. ((3j) , it is possible 
to rewrite the statistical model yicld in the coalcscence- 
likc form, 



, tat _ g h V H {m h T H f' 2 



N st at = 



3/T 



(2tí) 3 / 2 



11 9iV H (miT H )Vz 



(31) 



where we considcr the hadron h to be composcd of several 
constituents, B = J^. m* — rrih is the binding energy, 
and Ni t H represents the yield of the i-ih constituent at 
the volume Vh- This relation holds since the fugacity 
of a particle is equal to the product of the constituent 
fugacities, and the particle fugacity is related to the yield 
according to 



Ni 



H 



ghVH 



ih/T H í 



2n 



\ 3/2 

\m h T H ) 



(32) 



By using Eqs. ([21]) and ([15]). the ratio R% s is found to 



bc approximatcly given as 

M^ 2 {2T H /uj)^ n -^/ 2 {2T c /uj) L 



R CS =e -B/T H _ 



V C 



m 3 h /2 (í + 2T c /uj) n - 1 + 



-l+L 



— TT 



Ni/Ve 



Vh V N ilH /V H fjL (21, + 1) 



n 



(2/i)!l 



(33) 



Since Tq and Tu are the samé at small baryon chemical 
potcntial, the particle density Ní/Vc,Ní : h/Vh and the 
volume ratio Vc/Vh would be common at RHIC and 
LHC. As a result, the ratio of the production yield of 
hadrons in the coalescence model to that in the statistical 
model bccomes the samé at RHIC and LHC. 

In Fig. [2 we show the hadron mass dependence of the 
ratio R^ s in heavy ion collisions at RHIC. The results 
are similar in heavy ion collisions at LHC as discussed 
in the above. It is seen that the effect of the internal 
structures of exotic hadrons on R^ s is particularly large 
for light exotic hadrons. We also show the ratio for nor- 
mál hadrons. Here we háve included the resonance decay 
contribution to pseudoscalar mesons such as p — > 2tt in 
evaluating the ratio. This estimate would be closer to 
the observational condition. As already mentioned, the 
ratios for normál hadrons are found to be in the range of 
0.2-2. 
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It should be noted that the coalescence model may 
overcstimate the yield of very loosely bound molecules. 
Since the average hadron distance is around 2 fm and the 
temperature is Tp = 125 MeV at the thermal freeze-out, 
the loosely bound and spatially extended hadrons would 
dissociate easily through the final-state interactions with 
other hadrons. For example, the deuteron yield is calcu- 
lated to be around 1.4 per unit rapidity, which is largcr 
than the statistical model prediction (~ 0.3). Expcri- 
mental data seems to be consistent with the statistical 
model result, suggesting the possibility of later coales- 
cence of loose ly bound particles with a few MeV binding 
energies jl!8l |. 



VI. SUMMARY 



In this article, we háve proposed a new approach of 
studying cxotic hadrons in rclativistic hcavy ion colli- 
sions at RHIC and LHC. We háve considered the yiclds 
of proposed exotic hadrons; /o(980), ao(980), A'(1460), 
D s (2317), T^, AT(3872), Z+_(4430),_and T c ° 6 for exotic 
mesons, A(1405), 9+(1530), KKN, DN, D*N, 9 CS , BN 
and B*N for exotic baryons, H, KNN, ÍÍÍ2, H++, DNN 
and BNN for exotic dibaryons. To obtain the yields of 
these exotic multiquark hadrons or hadronic molecular 
states, we háve ušed the coalescence model based on ei- 
ther the quark degrecs of frccdom or the hadronic degrees 
of frecdom. 

Our results indicate that the yields of many exotic 
hadrons are large enough to be measurable in cxpcri- 
ments. In particular, heavy exotic hadrons containing 
charm and bottom quarks as well as strange quarks can 
be possibly observed at RHIC and especially at LHC. 
Therefore, relativistic heavy ion collisions will provide a 
good opportunity to search for exotic hadrons. In par- 
ticular, it may lead to the hrst observation of new ex- 
otic hadrons. Also, we háve found that the structure 
of light exotic hadrons has a signiheant effect on their 
yields in heavy ion collisions. For a hadron of normál 
quark structure, its production yield relative to the sta- 
tistical model prediction Rh = Nh/N^ at is found in the 
range of 0.2 — 2. The yield ratio is smaller (Rh < 0.2) if 
a hadron has a compact multiquark configuration. For a 
hadron of molecular configuration with an extended size, 
its yield is, on the other hand, largcr than the normál val- 
ues (Rh > 2). Therefore, the ratios of measured yiclds 
from experiments on heavy ion collisions to those pre- 
dicted by the statistical model provides a new method to 
diseriminate the different pictures for the structures of 
exotic hadrons. The study of exotic hadrons in relativis- 
tic heavy ion collisions thus will help answer longstand- 
ing problems on the existence and structure of exotic 
hadrons. 
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Appendix A: d-wave Wigner function 



In this Appendix, we extend the calculation shown in 
Ref. [l!9j | to construct the d-wave Wigner function from 
the harmonie oscillator wavc functions. For this, we need 
the wave function of the second excited statě, in addition 
to those of the ground and the first excited statě, as the 
basis functions for the Wigner function, 



Mx)= (é^y 






2 X _J=± 

— e 26^ 
b 



(j> 2 (x) 



ls)'0 



260F 



x 



1 e 26^ 



(AI) 



where b 2 = h/(mu>). The M-state wave functions of the 
L = 2 statě in Cartesian coordinates, are then 
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VffM?) 



\ý2(x)<l>o(y) - 0o(a;)02(y)j 



iV24>!(x)(j)i(y) 



<l>o(z) 



4>M=Ár) =^\<l>Áx)<fo{y) + i<fo{x)4>i{y))<i>\{z) 

- 4>2{x)4>o{y)(t>o{z) - 4>o{x)<t>2{y)<t>v{zyj 

<ftň-2<f) =<t> L uU^ (A2) 

From the dcnsity matrix of the L = 2 statě, which is 
defined by averaging over the M states, 

+ráf= 2 o (ř)<tirM?r + 0M=-i(r)^ri-i(^)* 

+^= 2 - 2 (r)0Íf = = 2 -2(O*), (A3) 

we obtain the M-averaged L = 2 harmonie oscillator 
Wigncr function, 



By replacing b with a(= 1/ \fJuJT), we obtain the d-wave 
coalescence factor using Eqs. (fTU| and (|A4j) as 



Jd 3 yd 2 kf(k)fW(y,k) 
fd 3 yd 2 kf(k) 



(47T ( T 2 )í 



2,2 



V(í + 2/xTcr 2 ) 3 15 



(2fj,Ta z ) 



(A5) 



Appendix B: Coalescence factor for the generál 
angular momentům l 

We derive in this Appendix the coalescence factor for a 
generál angular momentům ž-state given by Eq. (|17[) . For 
constituents that are uniformly distributed in space, we 
can integrate the Wigncr functions over space to obtain 
the momentům distribution of the constituents in the 
hadron. For the harmonie oscillator wave function of the 
lowest energy statě with a given /, it is given by 



P z (fc) = (W) 3 / 2 



(2a 2 fc 2 )< e _ a , k . 



(21 + IV! 

with the normalization J Pi(k)d 3 k/(2ir) 3 
lescence factor is then given by 



(Bl) 
1. The coa- 
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